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Abstract 

Evidence for the energy threshold of creating the quark-gluon plasma in nucleus- 
nucleus collisions, the so-called onset of deconfinement, has been found by the energy 
scan program of the NA49 experiment at the CERN SPS. In this paper we review 
the experimental and theoretical status of this phenomenon. First, the basic, quali- 
tative ideas are presented for non-experts. Next, the latest experimental results are 
compared to a statistical model within which the onset of deconfinement and its 
signals had been predicted. Finally, alternative interpretations and open questions 
are discussed. 
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1 Introduction 



One of the important issues of contemporary physics is the understanding of strong in- 
teractions and in particular the study of the properties of strongly interacting matter in 
equilibrium. What are the phases of this matter and what do the transitions between 
them look like ? These questions motivate broad experimental and theoretical efforts 
since more than 40 years. The study of high energy collisions between two atomic nuclei 
give us the unique possibility to address these issues in well controlled laboratory experi- 
ments. In particular, the advent of the quark model of hadrons and the development of 
the commonly accepted theory of strong interactions, quantum chromodynamics (QCD), 
naturally led to expectations that matter at very high densities may exist in a state of 
quasi-free quarks and gluons, the quark-gluon plasma (QGP) [Il[2l|3|. 

Experimental searches for QGP signals started at the Super Proton Synchrotron (SPS) 
of the European Organization for Nuclear Research (CERN) and the Alternating Gradient 
Synchrotron (AGS) of Brookhaven National Laboratory (BNL) in the mid 1980s. Today 
they are pursued also at much higher collision energies at the Relativistic Heavy Ion 
Collider (RHIC) at BNL. Soon experiments on nucleus-nucleus collisions at the Large 
Hadron Collider (LHC) in CERN will join the world effort at energies 20 times higher 
than at RHIC. Most probably that the QGP is formed at the early stage of heavy ion 
collisions at the top SPS energy and at RHIC energies. Unambiguous evidence of the 
QGP state was however missing. This may be attributed to the difficulty of obtaining 
unique and quantitative predictions of the expected QGP signals from the theory of strong 
interactions. 

For this reason the NA49 Collaboration at the CERN SPS has searched over the past 
years for signs of the onset of QGP creation in the energy dependence of hadron production 
properties. This search was motivated by a statistical model |1] showing that the onset 
of deconfinement should lead to rapid changes of the energy dependence of numerous 
experimentally detectable properties of the collisions, all appearing in a common energy 
domain. The predicted features have recently been observed ^ and dedicated experiments 
now continue detailed studies in the energy region of the onset of deconfinement. 

It is thus time for a summary. In this paper we review the experimental and theoretical 
status of the onset of deconfinement. First, the basic qualitative ideas are presented for 
non-experts. Next, a quantitative model within which the onset of deconfinement and its 
signals were predicted is reexamined and compared with the latest experimental results. 
Finally, alternative interpretations and open questions are discussed. 

2 Onset of Deconfinement for Pedestrians 

Phase transitions are fascinating physical phenomena. Small changes in temperature or 
pressure lead to dramatic changes in macroscopic properties of matter. Common examples 
from our daily life are transitions between solids, liquids and gases like boiling and freezing 
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of water. The well known phase diagram of water is shown in Fig. [T} where the regions 




Temperature (K) 

Figure 1: Phases of water. When adding heat (energy) at constant pressure water is 
transformed from solid to liquid and then from liquid to vapor as indicated by the dashed 
arrow. 

of existence of the various phases of water are depicted in a diagram of pressure and 
temperature. When adding heat to water one increases its temperature moving through its 
different phases and crossing their boundaries, as indicated by the dashed arrow in Fig. [T] 
for the example of constant atmospheric pressure. Dependence of the water temperature 
on the amount of added heat, called the heating curve of water, is shown in Fig. [2j In 
pure phases, such as ice, water or vapor, the temperature increases monotonically with 
added heat. The two regions of constant temperature (steps) signal the ice- water and 
water-vapor phase transitions. In these mixed phase regions added heat is used for the 
phase transformation instead of the increase of temperature as in the pure phase regions. 

Properties of water and other substances surrounding us and the transitions be- 
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Figure 2: Heating curve of water at fixed atmospheric pressure. It corresponds to the 
trajectory in the phase diagram of water indicated by the dashed arrow in Fig. [T] 



tween their various phases are determined by electromagnetic interactions of atoms and 
molecules. On the other hand, the properties of atomic nuclei which are built from nucle- 
ons (protons and neutrons) are determined by strong interactions. Naturally, the question 
arises whether strongly interacting matter also exists in distinct phases. What are their 
properties? At which temperatures do the transitions between them take place? What 
do these transitions look like? 

Since more than 40 years it is known that hadrons (i.e. mesons and baryons; all 
strongly interacting particles observed in nature are called hadrons) consist of more ele- 
mentary particles, the quarks and gluons. However, isolated quarks or gluons were never 
observed. They seem to be always confined in the interior of hadrons. But could a different 
phase of strongly interacting matter exist in which quarks and gluons are deconfined ? 

There are 3 parameters which describe the thermodynamical properties of a system. In 
non-relativistic systems they are temperature, particle number density, and pressure. The 
equation of state connects them, e.g., the pressure is a well defined function of temperature 
and particle density for a specific substance. In experiments on water one can most easily 
fix temperature and pressure to define the point on the phase diagram in Fig. [TJ Unlike in 
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water, the number of particles is not conserved in strongly interacting relativistic matter. 
Instead of particle number density the baryonic number, i.e. the difference between the 
number of baryons and anti-baryons, is conserved. In calculations it is convenient to use 
the equivalent variables baryonic number density or baryonic chemical potential. The 
phase diagram of strongly interacting matter emerging from theoretical considerations 
and experimental results is shown in Fig. |3] in terms of the commonly used variables 
temperature and baryonic chemical potential. Laboratory experiments (see discussion 
below) can create strongly interacting matter with different temperatures T and baryonic 
chemical potentials /i^. The functional dependence of the pressure on T and /i^, i.e. 
the equation of state of strongly interacting matter, remains the subject of intensive 
experimental and theoretical studies. 

A transition between the deconfined and the confined phase of strongly interacting 
matter probably took place during the expansion and cooling of the early Universe, about 
1 microsecond after the Big Bang. Cosmological signatures of this transition are difficult 
to identify today. However, extremely dense strongly interacting matter fills the interior 
of neutron stars. Arguments in favor of the existence of quark matter in the center of such 
stars were advanced already in the 1960s, soon after formulation of the quark hypothesis. 
One of the pioneering papers [2] argued: "A neutron has a radius of about 0.5-1 fm (1 fm 
= 10~^^ m), and so has a density of about 8-10^^ g-cm~^, whereas the central density of 
a neutron star can be as much as 10^^ — 10^^ g-cm~^. In this case, one must expect the 
hadrons to overlap, and their individuality to be confused. Therefore, we suggest that 
matter at such densities is a quark soup." The creation of matter in a deconfined phase, 
i.e. in the QGP phase, may be the only possibility to 'see' quarks and gluons moving 
freely in a large volume. 

Cosmological and astrophysical objects with the required properties are, unfortunately, 
difficult to investigate. Systematic study of the properties of strongly interacting matter 
requires a method to create it under well controlled conditions in the laboratory. The 
study of collisions of two heavy nuclei gives us this possibility. Such a collision produces 
a droplet of strongly interacting matter of high energy density, the so-called fireball. It 
is natural to expect that with increasing collision energy the fireball energy density also 
increases. Thus like in the case of water heating and observing successive transitions 
between its phases, we hope that with increasing collision energy we can detect anomalies 
in the energy dependence of hadron production properties and thus discover successive 
transitions between various phases of strongly interacting matter created at the early stage 
of collisions. The arrow in Fig. |3] schematically traces the position of the initially created 
fireball on the phase diagram when the energy of nucleus-nucleus collisions is increasing. 
At sufficiently high collision energy this matter droplet may reach the QGP phase (see 
Fig. [s]). Unfortunately the life time of the fireball is very short, about 10~^^ seconds. It 
quickly expands, cools down (see Fig. |4]) and finally decays into hadrons and a few light 
nuclei. These decay products are measured in detectors surrounding the collision point. 

The first phase transition of strongly interacting matter was observed studyingg col- 
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Figure 3: Phases of strongly interacting matter. With increasing coUision energy the 
matter created at the early stage of nucleus-nucleus collisions changes its properties as 
indicated by the arrow. At low energy it is in the confined phase (hadrons), at sufficiently 
high energy in the deconfined phase (QGP). M is the critical point of the nuclear liquid- 
gas phase transition. The shaded band shows the 1** order phase boundary between the 
hadron and QGP phase which is expected to end in a critical end point E. At E the 
sharp phase transition turns into a rapid crossover indicated as the dotted line. 



lisions at very low energies [6] (the energy per nucleon-nucleon pair in the center of mass 
system y/s^j^ < 2 GeV). This transition between a nuclear liquid and a nuclear gas hap- 
pens at a temperature of about 6T0^° K (5 MeV). The phase transition line and critical 
point M lie at large and small T inside the hadron phase region of the phase diagram 
as shown in Figs. [3] and |4j 

Emerging results from the study of high energy collisions of nuclei confirm the existence 
of the second phase transition in strongly interacting matter which was suggested by QCD. 
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(MeV) 

Figure 4: Parameters of strongly interacting matter created at the early stage of nucleus- 
nucleus interactions are shown by the full circles for central Pb+Pb (Au+Au) collisions 
at the top AGS energy (-y/sjviv ~ ^-^ GeV), intermediate SPS energy (i/s^^v ~ '^■^ GeV) 
and top SPS energy (v^^^v ~ GeV). The created fireball expands and cools along 
trajectories indicated by solid lines and decouples at the freeze-out points (full squares, 
triangles and star). 



It is the so-called deconfinement phase transition. Let us briefly present how one looked 
for this phase transition, and explain the main results. 

The search for the phase of deconfined strongly interacting matter, the QGP, already 
has a long history. It received a boost from the first acceleration of oxygen and sulfur nuclei 
at the CERN SPS in 1986 i^/s^J^ ^ 20 GeV) and of lead nuclei in 1994 (^s^^ ^ 17 GeV). 
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Figure 5: Heating curves of strongly interacting matter. Hadron production properties 
(see text for details) are plotted as a function of collision energy for central Pb+Pb 
(Au+Au) collisions (upper set of points) and p+p interactions (lower set of points) [5] 



Measurements from an array of experiments indicated that the critical energy density was 
probably exceeded and matter with unusual properties appeared to be formed in the early 
stage of the collisions [7] . A key problem was the identification of experimental signatures 
of QGP creation. Several signatures of the formation of a transient QGP state during the 
early stage of the collision had been proposed in the past [SI IS]- However, the uniqueness 
of these signatures came under renewed scrutiny and they were found not to be specific 
for the creation of QGP (see Appendix 8.1 for details). 

In the mid 1990s a study of results from experiments at CERN and the AGS at 
BNL (maximum energy a/s^v^ ~ 5.5 GeV) raised pUl ITT] intriguing questions concerning 
the energy dependence of hadron production between top AGS and SPS energies. In 
response to these questions a statistical model of the early stage of the collision process was 
proposed |1] in which an equation of state with a 1^* order phase transition was assumed. 
In this model the onset of deconfinement led to the prediction of a non-monotonic collision 
energy dependence of several hadron production properties. In particular, the model 
predicted a sharp maximum in the ratio of multiplicities of strange hadrons (the hadrons 
which contain s and s quarks) to pions (the lightest hadron) at the beginning of the 
transition region, at about v^tvtv ~ "^-^ GeV. This prediction triggered an extension of 
the experimental program at the SPS, the energy scan program [12]. Within this program 
head-on (central) collisions of two lead nuclei (Pb+Pb) were registered at several lower 
SPS energies (-\/sjvAf ~ 6.3, 7.6, 8.7 and 12.3 GeV) by the NA49 experiment. Other heavy 
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ion experiments at the SPS (NA45, NA50, NA57 and NA60) participated in selected runs 
of this program [13] . Final results, obtained mainly by the NA49 collaboration, confirm 
the qualitative expectations and the quantitative predictions of the model: rapid changes 
in properties of hadron production occur within a narrow energy range, a/s^^ = 7- 
12 GeV t5j. 

The most dramatic effect is seen in the energy dependence of the ratio of total particle 
yields of kaons and pions, (K+)/(7r+), in central Pb+Pb collisions which is plotted in 
Fig. jSj/e/t. Following a fast threshold rise the ratio passes through a sharp maximum in 
the SPS range and then seems to settle to a plateau value at higher energies. Kaons are 
the lightest strange hadrons and due to approximate isospin symmetry the {K^) yield 
counts about half of the strange (anti-) quarks produced in the coUisions and contained 
in the reaction products (see Appendix 8.2 for details). Thus Fig. [s] /e/t demonstrates 
that the fraction of strangeness carrying particles in the produced matter passes through 
a sharp maximum at the SPS in nucleus-nucleus collisions. This feature is not observed 
for proton-proton reactions as shown by the open dots in Fig. [5] left. 

A second important result is the stationary value of the apparent temperature T of 
mesons in central Pb+Pb collisions at SPS energies as shown in Fig. [5]ri(7/i)f:. In the 
fireball picture the apparent temperature is related to the local random motion of the 
particles and their collective expansion velocity in the direction transverse to the collision 
axis. 

Presently the sharp maximum and the following plateau in the energy dependence 
of the (K^) / (ti^) ratio has only been reproduced by the statistical model of the early 
stage. In this model it reflects the decrease in the ratio of strange to non-strange number 
of degrees of freedom when deconfinement sets in. The stationary value of the apparent 
temperature of mesons was predicted [HI [151 [IB] as a consequence of the constant 
pressure and temperature at the early stage of nucleus-nucleus collisions in the SPS energy 
range due to the coexistence of hadronic and deconfined phases. 

These results serve as evidence that the deconfinement phase transition in 
Pb+Pb collisions starts in the SPS energy range. The exciting and rich physics 
which can be studied in nucleus-nucleus collisions at the CERN SPS energies motivates 
ongoing and future experimental programs at the CERN SPS [TTJ [18], BNL RHIC [T9] . 
FAIR SIS [20] and JINR NICA [21]. 

3 Statistical Model of the Early Stage (SMES) 

The experimental search for the onset of deconfinement performed by the experiment 
NA49 at the CERN SPS was motivated by the predictions of the Statistical Model of the 
Early Stage [1], which treats the creation of the fireball in nucleus-nucleus collisions in a 
statistical model approach. The model does not attempt a description of the subsequent 
(hydro-) dynamical evolution of the fireball. First, this section reviews the main assump- 
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tions of the originally formulated model and presents results obtained from analytical and 
numerical calculations as well as comparison to the experimental data available 12 years 
ago, when the model was formulated. Next, later extensions of the SMES are discussed 
which address collective flow at freeze-out [I6l |22] and event-by-event fluctuations [231 121] • 
The subsequent Sec. 4 shows the comparison of the predictions of the model with the most 
recent experimental results. 

3.1 Main Assumptions 

1. The basic assumption of the SMES is that the production of new degrees of freedom 
during the early stage of A+A collisions is a statistical process. Thus formation of all 
microscopic states allowed by conservation laws is equally probable. As particle creation 
from energy does not produce net charges, only states with total baryon, flavor and electric 
charge quantum numbers equal to zero are considered. Presence of the colliding nucleons 
is assumed to affect the properties of the observed final state only via their interactions 
with the statistically produced particles during the expansion of the system. This issue 
is further discussed in point 5 below. Consequently, the properties of the state produced 
at the early stage are entirely defined by the available energy and the volume in which 
production takes place. In central A+A collisions this volume is chosen as the Lorentz 
contracted volume occupied by the colliding nucleons (participant nucleons) from a single 
nucleus: 



where Vq = ^tivqAp and 7 = y/sj^j^/ {2m]y), miy is the nucleon mass and Ap is the number 
of participant nucleons from a single nucleus. The tq parameter is taken to be 1.30 fm in 
order to fit the mean baryon density in the nucleus, po = 0.11 fm~^. 

2. Only a fraction, rj, of the total energy in A+A collision is transformed into the energy 
of new degrees of freedom created in the early stage. This is because a part of the energy 
is carried by the net baryon number which is conserved during the collision. The released 
(inelastic) energy can be expressed as: 

E = r]{^/s^^ -2mN) Ap . (2) 

The parameter r] is assumed to be independent of the collision energy and the system size 
for A+A collisions. The value of t] used for the numerical calculations is 0.67 [25] . 

3. The elementary particles of strong interactions are quarks and gluons. The deconfined 
state is considered to be composed of u, d and s quarks and the corresponding anti- 
quarks each with internal number of degrees of freedom equal to 6 (3 color states x 2 
spin states). The contribution of c, b and t quarks can be neglected due to their large 
masses. The internal number of degrees of freedom for gluons is 16 (8 color states x 2 
spin states). The masses of gluons and non-strange (anti-)quarks are taken to be 0, the 
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strange (anti-)quark mass is taken to be 175 MeV |26]. The properties of equilibrated 
matter is characterized by an equation of state (EoS). For the case of colored quarks and 
gluons the model assumes the ideal gas EoS modified by a bag constant B (see, e.g., 

mm)- 

p = p'''-B , e = e"' + B , (3) 

where p and e denote pressure and energy density, respectively, and the superscript ^'^ 
marks the quantities for the ideal gas. This equilibrium state is called the Quark Gluon 
Plasma or Q-state. 

4. The model uses an effective parametrization of the confined state, denoted as W- 
state (White-state). The non-strange degrees of freedom which dominate the entropy 
production are taken to be massless bosons, as suggested by the original analysis of 
entropy production in N+N and A+A collisions f28J. Their internal number of degrees 
of freedom was fitted to the data [28] and appeared to be about 3 times lower than 
the internal number of effective degrees of freedom in the QGP. The internal number of 
degrees of freedom for a QGP is 16 + (7/8) -36 = 48 and therefore the internal number 
of non-strange degrees of freedom for low energy collisions is taken to be 48/3 = 16. The 
mass of strange degrees of freedom is assumed to be 500 MeV, equal to the kaon mass. 
The internal number of strange degrees of freedom is estimated to be 14 as suggested by 
the fit to the strangeness and pion data at the AGS. Also for the W-state the ideal gas 
EoS is selected. Clearly, this description of the confined state should only be treated as an 
effective parametrization. The numerical parameters are fixed by fitting A+A data at the 
AGS and the parametrization is then used for extrapolation to higher collision energies 
where the transition between the confined and deconfined state is expected. 

5. It is assumed that the matter created at the early stage expands, hadronizes and 
freezes-out. Within the original SMES formulation these later stage were not modelled. 
It was, however, postulated that during these stages the total entropy and total number of 
s and s quarks created in the early stage are conserved. The only process which changes 
the entropy content of the produced matter during the expansion is assumed to be the 
interaction with the baryonic subsystem. It was argued that this leads to an entropy 
transfer to baryons which corresponds to the effective absorption of about 0.35 vr-mesons 
per baryon [29] . Thus the final hadronic state has non-zero baryonic number and electric 
charge. 

3.2 Analytical Formulas 

In this section the simplified version of the model (massless particles) will be discussed 
which allows to perform calculations analytically. Subsequently in Sec. 4 quantitative 
results from numerical calculations using finite masses will be presented and compared to 
measured data. All chemical potentials have to be equal to zero, as only systems with 
all conserved charges equal to zero are considered. Thus the temperature T remains the 
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only independent thermodynamical variable. It is convenient to define the EoS in terms 
of the pressure function p = p{T) as the entropy and energy densities can be calculated 
from the thermodynamical relations: 

In the case of an ideal gas the pressure of the particle species 'j' is given by: 

1 

(5) 



P'iT) = ^ / k'dk 



27r2 Jo 3(F + m2)V2 



exp I ^ I ± 1 



where is the internal number of degrees of freedom (degeneracy factor) for the j- 
th species, rrij is the mass of the particle, appears in Eq. ([s]) for bosons and '+1' 
for fermions. The pressure p{T) for an ideal gas of several particle species is additive: 
p{T) = 'Y^jP^iT). The same is valid for the entropy and energy densities of Eq. (|4). 

In order to be able to perform analytical calculations of the system entropy and illus- 
trate the model properties it is assumed that all degrees of freedom are massless. In this 
simplified case the pressure Eq. (|5| is equal to: 



pi{T) = jT\ (6) 

where is the so called Stephan-Boltzmann constant, equal to tt^q^ /30 for bosons and 
Ivr^fyf-'/SO for fermions. The total pressure in the ideal gas of several massless species can 
then be written as p{T) = n'^gT^ /90 with the effective number of degrees of freedom g 
given by 

9 = 9'+l9^, (7) 

where g'' and g-^ are internal degrees of freedom of all bosons and fermions, respectively. 
The g parameter is taken to be gw for the W-state and gQ for the Q-state, with gg > gy/. 
The pressure, energy and entropy densities then follow as: 

Pw(r) = %^rS .„,.(T) = ^1*: T' , .,,(T) = , (8) 

pq(t) T' - B , .g(r) = ^T« + B ,.g(r) = ^T\ (9) 

for the pure W- and Q-state, respectively. Note the presence of the non-perturbative bag 
terms in addition to the ideal quark-gluon gas expressions for the pressure and energy 
density of the Q-state. 
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The 1^* order phase transition between W- and Q-state is defined by the Gibbs crite- 
rion 

Pwin) = pq{T,) , (10) 
from which the phase transition temperature can be calculated as: 

905 



1/4 

(11) 



At T = Tc the system is in the mixed phase with the energy and entropy densities given 
by 

^mix = (1 ~ ^)^W ~^ ^^Q ' ^mix = (1 ^ ^)^W ~^ ^^Q ' i^'^) 

where (1 — ^ the relative volumes occupied by the W- and Q-state, respectively. 
From Eqs. (jsj [o]) one finds the energy density discontinuity ('latent heat') 

As = SQiT,) - swiT,) = e^Q-Ew = 45. (13) 

The early stage energy density is an increasing function of the collision energy and is 
given by (see Eqs. (Ill [2])): 



V 2mN ■ ^ ^ 

There is a remarkable equivalence (see Appendix C of Ref. [1]) of the Gibbs criterion 
(i.e. the pure phase corresponds to the larger pressure pw or pq and the mixed phase to 
equal pressures pw = Pq) and the maximum entropy criterion, 

s{e) = max { Swie), sq^e), SmiA^) } , (15) 

for an arbitrary EoS p = p{T) with a 1'^* order phase transition. For e < or e > Eq 
the system consists of pure W- or Q-state, respectively, with entropy density given by 
the following equations: 

M.) = ^(#V"(.-BP. (17) 



3 V 30 



For < e < Eq the system is in the mixed phase (Eq. 12) and its entropy density can 
be expressed as: 

Smixie) = — + e = a + be . (18) 



14 




100 



200 



300 400 
T [MeV] 



Figure 6: Energy density and pressure divided by T'^ as a function of temperature T. 
The bag constant B was adjusted to 600 MeV/fm^ to obtain a critical temperature 



Tc = 200 MeV. 



The ratio of the total entropy of the created state to the number of nucleons participating 
in A+A collisions is: 

S V s rriN s 



2Ar 



2A 



and is independent of the number of participant nucleons. The entropy density s in 



Eq. (19) is given by the general expressions Eq. (15) with e defined by Eq. (14). For small 
y/sj^j^ the energy density Eq. ( 14 ) corresponds to that of the pure W-state and one finds 



S 



2Ar 



where 



C 



2 f Tx'^mi^\^^^ 
15po 



p/ w 



(20) 



(21) 
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Thus for low collision energies, where the W-state is created, the entropy per participant 
nucleon is proportional to F. For high a/sjvat pure Q-state is formed and Eq. (19) 
leads to 

5" \ ^ i/i ( 2mNB 



Cg'A-F [I - ^ ^ r. (22) 



C F ( 1 '"^^^ 



2r]poF^ 



For large values of F the entropy per participant nucleon in the Q-state is also proportional 
to F. The slope is, however, larger than the corresponding slope for the W-state by a 
factor (gq/ gwY^'^- In the interval of F in which the mixed phase is formed the energy 
dependence of the entropy per participant nucleon is given by: 

(t^) = + (v^iviv - 2m^) , (23) 

\'^-^V/ mix V^NN 

where 

Ci = ^ a , C2 = vb. (24) 
Po 



Equation (23) gives approximately a quadratic increase with F of the entropy per partic- 
ipant nucleon in the mixed phase region. 

Let us now turn to strangeness. The model defines g^ and gQ as the numbers of inter- 
nal degrees of freedom of (anti-)strangeness carriers in the W- and Q-state, respectively. 
The total entropy of the considered state is given by the sum of entropies of strange and 
non-strange degrees of freedom. Provided that all particles are massless the fraction of 
entropy carried by strange (and anti-strange) particles is proportional to the number of 
strangeness degrees of freedom: 

= ^ S . (25) 

g 

Equation (25) is valid for both W- and Q-state. Note that all degeneracy factors are 
calculated according to the general relation Eq. ([T]). For massless particles of the j-th 
species the entropy is proportional to the particle number: 

Sj = 4iV, . (26) 

Thus the number of strange and anti-strange particles can be expressed as 



iVs + iV^ = f - , (27) 

4 g 



and the strangeness to entropy ratio is equal to 

Ns + N^ ^ 1 / 
S 4 g 



(2^ 
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Figure 7: Left: The early stage (initial) temperature of the fireball as a function of F. 
Right: The fraction of volume occupied by the QGP as a function of F. 



One concludes therefore that the strangeness to entropy ratio for the ideal gas of massless 
particles is dependent only on the ratio of strange to all degrees of freedom, g^/g. This 
ratio is expected to be equal to gq/gq = 0.22 in the Q-state and gw/9w — 0.5 in the 
W-state. Consequently the phase transition from the W- to the Q-state should lead to 
a decrease of the strangeness to entropy ratio by a factor of about 2. This simple picture 
will be modified because of the large value of the mass of strange degrees of freedom in 
the W-state [m^ = 500 MeV) compared to the temperature T. In this case the left hand 
side of Eq. ( [28| is a strongly increasing function of T. 

In order to demonstrate properties of the EoS the ratios of e/T"^ and p/T^ are plotted 
in Fig. [6] as functions of the temperature. The bag constant B = 600 MeV/fm^ was 
adjusted such that the resulting critical temperature Tc is equal to 200 MeV. 

3.3 Quantitative Calculations 

We next turn to the results from numerical calculations based on the model using finite 
strangeness carrier masses. For the number of non-strange degrees of freedom gg^ and g^ 
one gets as in the simplified model: 

= 2 ■ 8 + ^ ■ 2 ■ 2 ■ 3 ■ 2 = 37 ; g"^ = IQ . (29) 
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The strange degrees of freedom are now considered to have reahstic masses rriQ and m^. 
Equation ^ is used with 

= 2 ■ 2 ■ 3 = 12 , ^ 175 MeV ; c/^ = 14 , ^ 500 MeV . (30) 

Note that there is no factor '7/8' in the above expression for gg as Eq. (Isj) with Fermi 
momentum distribution was taken. The contributions of strange degrees of freedom to the 
entropy and energy densities are calculated using the thermodynamical relations Eq. (|4]). 
A convenient variable to study collision energy dependence is the Fermi-Landau vari- 



able F defined in Eq. (21). The dependence of the early stage temperature T (initial 
temperature of the fireball) on F in the SMES is shown in Fig. [T] left. Outside the tran- 
sition region T increases approximately linearly with F. Inside the transition region T 
is constant (T = = 200 MeV). The transition region begins at F = 2.23 GeV^/^ 
(p^^^ = 30A GeV) and ends at F = 2.90 GeV^/^ ^p^^^ ^ Q^yy rj^^g fraction of the 
volume occupied by the Q-state, ^, increases rapidly in the transition region, as shown in 
Fig. [7] right. 
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Figure 8: Left: The entropy per participant nucleon as a function of F (solid line). 
Dashed line indicates the dependence obtained assuming that there is no transition to the 
QGP. Right: Ratio of produced entropy in pion units per participant nucleon, {St,)/ (Np), 
as a function of F. Experimental data on central collisions of two identical nuclei are 
indicated by closed circles. These data correspond to the status of 1998 [TOl [11] and 
should be compared with the model predictions shown by the solid line. The open boxes 
show results obtained for nucleon-nucleon interactions. 
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The number of non-strange and strange degrees of freedom and their masses are given 



by Eqs. (29, 30). They are different in the confined and deconfined phases. Thus, one 



expects abrupt changes of the pion multiphcity (entropy) (see Fig. [8] left) and the multi- 
phcity of strange particles (see Fig. |9] left) as a function of colhsion energy in the energy 
range where a transition from confined to deconfined matter takes place at the early 
stage of A+A collisions. The comparison of these predictions with experimental results 
is discussed in Sees. 4.1 and 4.2. 




Figure 9: Left: The ratio of the total number of s and s quarks and anti-quarks to 
the entropy (solid line) as a function of F. The dashed line indicates the corresponding 
ratio calculated assuming absence of the phase transition to the QGP. Right: The ratio 
Es of strange particle to pion production as a function of F. Experimental data on 
central collisions of two identical nuclei are plotted as closed circles. These data should 
be compared with the model predictions shown by the solid line. The open boxes show 
results obtained for nucleon-nucleon interactions, scaled be a factor 3.6 to match A+A 
data at AGS energy. The plotted data show the status of 1998 as compiled in flUi ITlj . 



3.4 Extensions of the SMES 

Several extensions of the SMES were developed over the past 10 years which allow to 
predict signals of the onset of deconfinement related to the collective flow at freeze- 
out [m [22] and to event-by-event fluctuations [23], [2l]. These extensions are briefly 
presented in this subsection. 
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3.4.1 Collective Flow at Freeze-out 

The collective flow of matter at freeze-out depends on the properties of the early stage 
as well as on the expansion dynamics itself. Within the SMES the collision energy de- 
pendence of the early stage properties is predicted. In particular, in the energy range in 
which the mixed phase is created the pressure and temperature are constant and at the 
end of the mixed phase domain the pressure to energy density ratio reaches its minimum 
(the softest point of the EoS). From general hydrodynamic considerations this is expected 
to lead to a reduction of the buildup of transverse [TB] and longitudinal [22] collective flow 
at freeze-out. The corresponding signals are discussed in Sees. 4.3 and 4.4. 

3.4.2 Event-by- Event Fluctuations 

Up to this point only quantities averaged over many collisions (events) were considered. 
Next an extension of the SMES is reviewed which leads to predictions of fluctuations from 
event to event. 

The key additional assumption is that when the collision energy is fixed, the energy, 
which is used for particle production (inelastic energy) can still fluctuate. These dynamical 
energy fluctuations lead to dynamical fluctuations of macroscopic properties X of the 
matter, like its entropy and strangeness content [23]. The relation between them is given 
by the EoS. For example, different values of the energy of the early equilibrium state 
lead to different, but uniquely determined, entropies. Since the EoS shows an anomalous 
behavior in the phase transition region, this anomaly should also be visible in the ratio 
of entropy to energy fluctuations [23] . 

According to the first and the second principles of thermodynamics the entropy change 
6S is given as T6S = SE + p5V. For central A-l-A collisions, one expects 6V = 0. Within 
the SMES the ratio of entropy to energy fluctuations can then be calculated and expressed 
as a simple function of the p/e ratio [23] : 



Within the SMES model, confined matter (which is modelled as an ideal gas) is created at 
the early collision stage below a collision energy of 30^4 GeV. In this domain, the ratio p/e, 
and consequently the ratio, are approximately independent of the collision energy and 
equal about 1/3 and 0.56, respectively. The SMES model assumes that the deconfinement 
phase transition is of the first order. Thus, there is a mixed phase region, corresponding 
to the energy interval 30^-60 A GeV. At the end of this region the p/e ratio reaches a 
minimum (the 'softest point' of the EoS [IS]). Thus in the transition energy range the Re 
ratio increases and reaches its maximum, R^ ^ 0.8, at the end of the transition domain. 
Further on, in the pure deconfined phase, which is represented by an ideal quark-gluon 
gas under bag pressure, the p/e ratio increases and again approaches its asymptotic value 




(31) 
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1/3 at the highest SPS energy of 160A GeV. The numerically calculated predictions of 
the SMES are plotted in Fig. 10 left. The early stage energy and entropy fluctuations 
entering Eq. (31) are not directly observable, however, as argued in Ref. [23], they can 



be inferred from the experimentally accessible information on the final state energy and 
multiplicity fluctuations. 
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Figure 10: The collision energy dependence of the fluctuation signals of the onset of the 
deconfinement calculated within the SMES. Left: The shark fin in of the ratio of entropy 
to energy fluctuations Re (31) (see Ref. [23J). Right: The tooth structure in the ratio of 
strangeness and entropy fluctuations Rg/e (33) (see Ref. |2^). 



In Ref. [23] the energy dependence of dynamical strangeness fluctuations caused by 
dynamical energy fluctuations was studied within the SMES model. Defining A^^ as the 
total number of strange quark-anti-quark pairs created in an A+A collision one calculates 
the fluctuation ratio as: 

For T — )► oo the system is in the QGP phase. Strange (anti-)quarks can be considered 
as massless and the bag constant can be neglected. Then e oc and oc and 
consequently de/e = A-dT/T and drts/ng = 3-dT/T, which results in Rg = (3/4)^ = 0.56. 
In the confined phase, T < Tc, the energy density is still approximately proportional 
to due to the dominant contributions of non-strange hadron constituents. However, 
the dependence of the strangeness density on T is governed by the exponential factor. 
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Us oc exp(— m^), as T << ms = = 500 MeV. Therefore, at small T one finds 
de/e = 4-dT/T and dug/ns = rris-dT/T'^, so that the ratio Rg = rris/i'iT) decreases with 
T. The strangeness density is small and goes to zero at T — )■ 0, but the fluctuation 



ratio Rs Eq. (32) is large and increases to infinity in the zero temperature limit. One finds 
a non-monotonic energy dependence of Re with a maximum at the boundary between the 
mixed phase and the QGP [23]. A pronounced minimum-structure is expected in the 
dependence of Rg on the collision energy [21]. It is located at a/s^vat ~ 7-12 GeV 
{30 A — 60 A GeV), where the mixed phase is created at the early stage of A+A collision. 

Both entropy and strangeness fluctuation measures, Rg and Rg, show anomalous be- 
havior in the transition region: a maximum is expected for Re and a minimum for Rg. 
Consequently, an even stronger anomaly is predicted for the ratio: 

H - - rs3^ 



shown in Fig. 10 right. Experimental measurements of Rg/e may be easier than the 
measurements of Re and Rg because the ratio Rg/e requires measurements of particle 
multiplicities only, whereas both Re and Rg involve also measurements of particle energies. 

These predictions are discussed in Sec. 4.5 in the context of existing experimental 
data. 



4 Signals of the Onset of Deconfinement 

Next the predictions of the SMES model reviewed in Sec. 3 will be related to directly 
measurable qunatities and compared with available experimental results. In particular, 
their significance as evidence for the onset of deconfinement will be discussed in detail. 



4.1 The Kink 

The majority of all particles produced in high energy interactions are pions. Thus, pions 
carry basic information on the entropy created in the collisions. On the other hand, 
entropy production should depend on the form of matter present at the early stage of 
collisions. Deconfined matter is expected to lead to a final state with higher entropy than 
that created by confined matter. Consequently, it is natural to expect that the onset of 
creation of deconfined matter should be signaled by an enhancement of entropy and thus 
pion production. This simple intuitive argument can be quantified within the SMES. 

Equilibration between newly created matter and baryons is assumed to take place 
during the evolution of the system. It was argued that this equilibration causes transfer 
of entropy from the produced matter to baryons. The analysis of the pion suppression 
effect at low collision energies indicates that this transfer corresponds to the effective 
absorption of about 0.35 pion per participant nucleon [21]. It is further assumed that 
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there are no other processes which change the entropy content of the state produced 
during the early stage. 

For the comparison with the model it is convenient to define the quantity: 

{S^) = {n) + K{K + K)+a{Np) , (34) 

where (vr) is the measured total multiplicity of final state pions and {K + K) is the 
multiplicity of kaons and anti-kaons. The factor k = 1.6 is the approximate ratio between 
mean entropy carried by a single kaon to the corresponding pion entropy at chemical 
freeze-out. The term a{Np) with a = 0.35 is the correction for the discussed partial 
transfer of the entropy to baryons. The quantity (S"^) can thus be interpreted as the early 
stage entropy measured in pion entropy units. The conversion factor between S and {St,) 
is chosen to be 4 entropy units per pion at chemical freeze-out). 

The dependence of the entropy per participant nucleon on F is shown in Fig. |8] left. 
Outside the transition region the entropy increases approximately proportional to F, but 
the slope in the Q-state region is larger than the slope in the W-state region. 

The number of baryons which take part in the collision {2Ap in the model calculations) 
is identified now with the experimentally measured number of participant nucleons, (Np) . 
The fraction of energy carried by the produced particles is taken to be ?7 = 0.67 [25] and 
is assumed to be independent of the size of the colliding nuclei and the collision energy. 

The comparison made in 1998 between the ratio (Stt) / (Np) calculated from available 
measurements and the model is shown in Fig. |8] right. The parametrization of the W-state 
was chosen to fit the AGS data and, therefore, the agreement with low energy A+A data 
is not surprising. On the other hand the description of high energy (SPS) results obtained 
by the NA35 and NA49 Collaborations is essentially parameter free. 

The pion multiplicity is proportional to the initial entropy, and the (tt) / (Np) ratio 
can thus be calculated outside the transition region as, 

oc g'/' F , (35) 



{Np) 



where g = g^ = 16 for the initial state in the confined phase and g = gq = 47.5 for 
the initial state in the deconfined phase at T >> rrig. Therefore, the {tt)/{Np) ratio 
increases linearly with Fermi's energy measure F outside the transition region, and the 
slope parameter is proportional to g^^^ |28]. In the transition region, a steepening of 
the increase of pion production with energy is predicted, because of the activation of the 
partonic degrees of freedom. 

The recent compilation of data [5J on pion multiplicity produced in central Pb+Pb 



(Au+Au) collisions and p + p{p) interactions is shown in Fig. 11 which displays the mean 
pion multiplicity (vr) = 1.5 ((vr^) + (tt^)) per wounded nucleon as a function of F. The 
results from p + p{p) interactions are shown by the open symbols. Up to the top SPS 
energy the mean pion multiplicity in p + p interactions is approximately proportional 
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Figure 11: Energy dependence of the mean pion multiplicity per participant nucleon 
measured in central Pb+Pb and Au+Au collisions (full symbols), compared to the corre- 
sponding results from p + p{p) reactions (open symbols). The compilation of data is from 
Ref. [5]. 



to F. A fit of {7t)/{Np) = b ■ F yields a value of 6 ^ 1.063 GeV'^/^. For central 
Pb+Pb and Au+Au collisions (filled symbols in Fig. 11) the energy dependence is more 
complicated. Below 40^4 GeV {y/sj^j^ = 8.7 GeV) the ratio (n) / {Np) is lower in A+A 
collisions than in p+p{p) interactions (pion suppression) while at higher energies this ratio 
is larger in A+A collisions than in p + p{p) interactions (pion enhancement). A linear fit, 
{tt)/{Nw) = a + 6 ■ F for F < 1.85 GeV^/^ gj^^g ^ ^ _o.45 and b = 1.03 GeV~^/^. The 
slope parameter fitted in the range F > 3.5 GeV^^^ is 6 = 1.33. This is shown by the solid 
line in Fig. 11 (the lowest data point at the top RHIC energy was excluded from the fit). 
Thus, in the region 15y4-40y4 GeV between the highest AGS and the lowest SPS energy 
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the slope increases by a factor of about 1.3. This agrees with the SMES result: 




) 



1/4 




1/4 



1.31 , 



(36) 



where the increase is caused by the creation of a transient state of deconfined matter in 
the early stage of the collisions at energies higher than i/s^jy ~ '^■^ GeV (30^4 GeV). 



The enhanced production of strangeness was considered by many authors as a potential 
signal of QGP formation [8l|30l[3T]. The line of arguments is the following. One estimates 
that the strangeness equilibration time in the QGP is comparable to the duration of 
the collision process (< 10 fm/c) and about 10 times shorter than the corresponding 
equilibration time in hadron matter. It is further assumed that in the early stage of the 
fireball the strangeness density is much below the equilibrium density, e.g., it is given by 
the strangeness obtained from the superposition of nucleon-nucleon interactions. Thus it 
follows that during the expansion of the matter the strangeness content increases rapidly 
and approaches its equilibrium value provided matter is in the QGP state. In hadron 
matter the modification of the initial strangeness content is less significant due to the 
long equilibration time. This leads to the expectation that strangeness production should 
rapidly increase when the energy threshold for the production of deconfined matter is 
crossed from below. 

In the SMES the role of strangeness is different. This is because statistical production 
of particles is postulated and therefore also strange particles are assumed to be produced 
in equilibrium. Consequently possible secondary processes do not modify its value. At 
T = Tc the strangeness density is lower in the QGP than in confined matter. Thus, a 
suppression of strangeness production is expected to occur when crossing the transition 
energy range from below. The low level of strangeness production in N+N interactions 
as compared to the higher strangeness yield per participant nucleon in central A+A 
collisions (called strangeness enhancement) can be understood as mostly due to the effect 
of strict strangeness conservation (canonical suppression) imposed on the strange and 
anti-strange degrees of freedom [32]. This constraint has an important effect for small 
statistical systems such as the confined matter in the early stage of N+N collisions. 

We are interested in the collision energy region between the AGS and SPS. At 'low' 
collision energies (when a pure W-state is formed) the strangeness to entropy ratio in- 
creases with F. This is due to the fact that the mass of the strange degrees of freedom 
is significantly higher than the system temperature. At T = Tc the ratio is higher in the 
W-state than in the Q-state region. Therefore the ratio decreases in the mixed phase 
region to the level characteristic for the Q-state. In the Q-state, due to the low mass of 
strange quarks compared to the system temperature, only a weak dependence of the ratio 



4.2 The Horn 
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on F is observed. The F dependence of the strangeness/entropy ratio calculated in the 
SMES is shown in Fig. |9] left. 

The comparison between the model and the data on strangeness production is per- 
formed under the assumption that the strangeness content created during the early stage 
is preserved till the hadronic freeze-out. This simplifies the model calculations by ne- 
glecting possible gluon contribution to strangeness production during hadronization of 
the QGP. 

Total strangeness production can be studied (in the AGS and SPS energy range) using 
the experimental ratio: 

E. ^ M±Z) . (3.) 

(tt) 

where (A) is the mean multiplicity of A hyperons (see Appendix 8.2). Within the SMES 



model Es of Eq. (37) is calculated as: 



F - ^^s + N,)/C 

- (5-5.)/4-a(Ar,) ' ^^S) 

where C, = 1.36 is the experimentally estimated ratio between total strangeness production 
and strangeness carried by A hyperons and K + K mesons [33] and Sg is the fraction of the 
entropy carried by the strangeness carriers. The comparison between the calculations and 
the data available in 1998 is shown in Fig.[9]n5'/ii [1]. The good description of the AGS data 
is again a consequence of the parametrization of the W-state: g^r = 14, m^^ = 500 MeV 
which was based on these data. As in the case of the pion multiplicity, the description 
of the strangeness results at the SPS (from the NA35 and NA49 Collaborations) can 
be considered as being essentially parameter free The agreement with the SPS data 
is obtained assuming creation of globally equilibrated QGP in the early stage of A+A 
collisions. The characteristic non-monotonic energy dependence of the Es ratio is a 
signature of the phase transition. 

The entropy and strangeness production in central A+A collisions satisfies well the 
conditions needed for thermodynamical treatment. Therefore one expects that the mea- 
sures of the entropy per participant nucleon, (S'^)/(iVp), and the ratio of strangeness to 
entropy production, Es, are independent of the number of participants for large enough 
values of {Np). 

The energy dependence of the strangeness to entropy production ratio is a crucial 
signal of deconfinement. The temperature dependence of the multiplicity of a particle is 

^ The Es value resulting from a QGP can be estimated in a simple way. Assuming that = 0, and 



neglecting the small (< 5%) effect of pion absorption at the SPS, one gets from Eq. (28) and Eq. (38) 
Es ~ ((?Q/1.36)/gQ* « 0.21, where Qq ~ (7/8) • 12 is the effective number of degrees of freedom of s and 
s quarks and gg' = 16 + (7/8) • 24 is the corresponding number for u, u, d, d quarks and gluons. Moreover, 
we use the approximation that the pion entropy at freeze-out is equal to the mean entropy of q, q and g 
in a QGP. 
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strongly dependent on its mass. In the Boltzmann approximation one finds: 



W = ^ Jpdpexpl 1 = —m^TK2[—) , (39) 

^ ^ 

where K2 is the modified Hankel function. For hght particles {mi/T ^ 1) one finds 



from Eq. (39), {Ni) oc T^, whereas for heavy particles {rrih/T ^ 1) Eq. (39) leads to 
{Nh) oc T^exp{— rrih/T). Within the SMES the strangeness to entropy production 
ratio increases steeply at low collision energies, when confined matter is produced. This 
is due to the low temperature at the early stage and the high mass of the carriers of 
strangeness (the kaon mass). Thus, uik ^ T and total strangeness production is pro- 
portional to T^/^exp(— rriK/T). On the other hand, the total entropy is approximately 
proportional to T^. Therefore, the strangeness to pion production ratio is approximately 
T~^/^ exp(— mx/T) in the confined phase and strongly increases with the collision energy. 
When the transition to deconfined matter occurs, the mass of the strangeness carriers is 
significantly reduced (m^ = 175 MeV, the strange quark mass). Due to the lower mass 
{rris < T) the strangeness yield becomes approximately proportional to the entropy (both 
are proportional to T^), and the strangeness to entropy (or pion) production ratio becomes 
independent of energy in the QGP. This leads to a 'jump' in the energy dependence from 
the larger value for confined matter to the value for deconfined matter. Thus, within the 
SMES, the non-monotonic energy dependence of the strangeness to entropy production 
ratio followed by a plateau at the deconfined value is a direct consequence of the onset of 
deconfinement taking place at y/sj^^ ~ 7.6 GeV (about 30 A GeV) [1]. 

The Es ratio was the first observable used to establish the energy dependence of the 
strangeness to entropy production ratio from the data. This ratio is closely proportional 
(see Appendix 8.2) to the {K^)/{7r^) ratio, which with time became better measured 



experimentally. The energy dependence of both ratios is plotted in Fig. [12] for central 
Pb+Pb (Au+Au) collisions and p + p interactions as function of collision energy. 

For p + p interactions the ratios show a monotonic increase with energy. In contrast, 
a very different behavior is observed for central Pb+Pb (Au+Au) collisions. The steep 
threshold rise of the ratios characteristic for confined matter changes at high energy into 
a constant value at the level expected for deconfined matter. In the transition region 
(at low SPS energies) a sharp maximum is observed caused by the higher strangeness to 



entropy production ratio in confined matter than in deconfined matter. As seen in Fig. 12 
the measured dependence is consistent with that predicted within the SMES 



4.3 The Step 

The energy density at the early stage increases with increasing collision energy. At low 
and high energies, when a pure confined or deconfined phase is produced, this leads to 
an increase of the initial temperature and pressure. This, in turn, results in an increase 
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Figure 12: Left: Energy dependence of the {K^) / {-k^) ratio measured in central Pb+Pb 
and Au+Au collisions (upper set of symbols) compared to the corresponding results from 
p+p reactions (lower set of symbols). Right: Energy dependence of the relative strangeness 
production as measured by the Es ratio (see text) in central Pb+Pb and Au+Au collisions 
(upper set of symbols) compared to results from p + p reactions (lower set of symbols). 
The compilation of data is from Ref. [5J. The dashed-dotted line in the figure shows the 
predictions of the SMES 



of the transverse expansion of the produced matter and consequently a flattening of the 
transverse mass spectra of final state hadrons. One may expect an 'anomaly' [HI UHl US] 
in the energy dependence of the average hadron transverse momenta in the mixed phase 
region where the temperature and pressure are approximately constant. 

The experimental data on spectra of the transverse mass rriT = {rn^+p^Y^'^ are usually 
parameterized by a simple exponential dependence: 



dN 



iTLTdrnT 



C exp (- . (40) 



The inverse slope parameter T* is sensitive to both the thermal and collective motion 
in the transverse direction. Hydrodynamical transverse flow with collective velocity vt 
modifies the Boltzmann mr-spectrum of hadrons. At low transverse momenta, it leads 
to the result {Tkin is the kinetic freeze-out temperature): 



^/L-p, = Tui^ + \mvl. (41) 
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Such a linear mass dependence of T* is supported by the data for hadron spectra at 
small pt- However, for pt ^ rn the hydrodynamical transverse flow leads to the mass- 
independent blue-shifted 'temperature': 



Thigh-pT - '^kin ■ \l ■ (42) 



Note that a simple exponential fit Eq. (40) neither works for light vr- mesons, T^ow-pt(^) ^ 
Thigh-pT('^)y nor for heavy (anti-)protons and ( ant i-) lamb das, T;;^_p^(p,A) > 
TLgH-pM^) (see e.g., Refs. [SIES]). 

Kaons are the best suited among measured hadron species for observing the effect 
of the modification of the EoS due to the onset of deconfinement in hadron transverse 
momentum spectra. The arguments are the following. First, the kaon mr-spectra are 
only weakly affected by hadron re-scattering and resonance decays during the post- 
hydrodynamic hadron cascade at SPS and RHIC energies jM]. Second, a simple one 



parameter exponential fit Eq. (40) is quite accurate for kaons in central A+A collisions 
at all energies. This simplifies the analysis of the experimental data. Third, high quality 
data on m-r-spectra of and K~ mesons in central Pb+Pb (Au-|-Au) collisions are 
available over the full range of relevant energies. 
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Figure 13: Energy dependence of the inverse slope parameter T* of the transverse mass 
spectra of K'^ (left) and K~ mesons (right) measured at mid-rapidity in central Pb+Pb 
and Au+Au collisions. The slope parameters are compared to those from p + p 
reactions (open circles). The compilation of data is from Ref. [3]. 
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Figure 14: Energy dependence of the mean transverse mass, (mr) — m, measured at 
mid-rapidity in central Pb+Pb and Au+Au collisions for tt^ {left), {middle), and p 
and p {right). Results for positively (negatively) charged hadrons are shown by full (open) 
symbols. The compilation of data is from Ref. 0. 



The experimental results on the energy dependence of the inverse slope parameter of 
and K~ transverse mass spectra for central Pb+Pb (Au+Au) collisions are shown 



in Fig. 13 The striking features of the data can be summarized and interpreted [16] as 
follows. The T* parameter increases strongly with collision energy up to the SPS energy 
point at y/sj^j^ = 7.6 GeV (30^4 GeV). This is the energy region where the creation 
of confined matter at the early stage of the collisions is expected. Increasing collision 
energy leads to an increase of the early stage temperature and pressure. Consequently 
the transverse momenta of produced hadrons, measured by the inverse slope parameter, 
increase with collision energy. This rise is followed by a region of approximately constant 
value of the T* parameter in the SPS energy range y/sj^j^ = 7.6 - 17.2 GeV {30A - 
158 A GeV), where one expects the transition between confined and deconfined matter 
with the creation of mixed phases. The resulting modification of the equation of state 
'suppresses' the hydrodynamical transverse expansion and leads to the observed plateau 
structure in the energy dependence of the T* parameter [16]. At higher energies (RHIC 
data), T* again increases with the collision energy. The EoS at the early stage becomes 
again stiff and the early stage pressure increases with collision energy, resulting in a 
resumed increase of T*. As also shown in Fig, 13 the parameter T* in p + p interactions 
appears to increase smoothly and does not show the step-like structure. 

For the transverse mass spectra of pions and protons the inverse slope parameter de- 
pends on the transverse mass interval used in the fit. The mean transverse mass {rrix) 
provides an alternative characterization of the my-spectra. The energy dependence of 
(m^) — m for pions, kaons and (anti-)protons is shown in Fig, 14 These results demon- 
strate that the approximate energy independence of (mr) — m in the SPS energy range 
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is a common feature for all abundantly produced particle species. 



4.4 The Dale 

This subsection summarizes the analysis of Ref. [22] of the longitudinal pion spectra within 
Landau's hydrodynamical model [36l |37]. The interest in this model was revived by the 
remarkable observation that the rapidity distributions at all investigated energies can be 
well described by a single Gaussian (see [3H] and references therein). Moreover, the energy 
dependence of the width can also be described reasonably well by the same model. 

The main physics assumptions of Landau's picture are as follows. The collision of 
two Lorentz-contracted nuclei leads to a complete stopping of the colliding nuclei and full 
thermalization of the created hadronic matter. This establishes the volume and energy 
density for the initial conditions of hydrodynamic expansion at each collision energy. 
Assuming for simplicity the equation of state in the form p = c^e {cg denotes the speed 
of sound, and = 1/3 for an ideal massless particle gas) the pion rapidity spectrum is 
given by |39lll0]: 

with 



,2 



cr, 



2 » 



y 31-cf 



ln(v^NN/2mN) , (44) 



where i^' is a normalization factor converting entropy to pion density The above pre- 
diction was compared with the experimental data on rapidity distributions of negatively 
charged pions produced in central Pb+Pb (Au+Au) collisions at various energies. Fig- 
ure 15 left shows the measured width cTy of the rapidity spectra [HI ESI 1121 US] as a function 
of the beam rapidity. The full line shows a linear fit through the data points. The dot- 
ted line indicates the Landau model predictions with = 1/3. The model only roughly 
reproduces the measured dependence. At low AGS energies and at the top RHIC energy, 
the experimental points are under-predicted, while in the SPS energy regime Landau's 
model over-predicts the width of the rapidity distributions. 

These deviations can be attributed to the changes in the EoS, which can be effectively 
parameterized by allowing the speed of sound to be dependent on collision energy. By 



inverting Eq. (44) one can express cj. in the medium as a function of the measured width 



^ There are two issues related to derivation of Eq. (44) which need clarification by future study. First, 



Eq. ( 44 ) is obtained assuming that Cg depends only on the early stage energy density and its dependence 
on decreasing energy density during expansion is neglected. Second, the Landau model assumes stopping 
and thermalization of the total energy in the cm. system, whereas only a fraction of the inelastic energy 
is stopped and thermalized in the SMES model. 
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Figure 15: Comparison of the Landau hydro dynamic model with rapidity distributions of 
charged particles [22]. Left: The root mean square width ay of the rapidity distributions 
of negatively charged pions in central Pb+Pb (Au+Au) reactions as a function of the 
beam rapidity i/p. The dotted line indicates the Landau model prediction with cj. = 1/3, 
while the full line shows a linear fit through the data points. Data (full symbols) are 
taken from [381 H2[ H3] . The statistical errors given by the experiments are smaller 
than the symbol sizes. Systematic errors are not available. Right: Speed of sound as 
a function of beam energy for central Pb+Pb (Au+Au) reactions as extracted from the 
data using Eq. (45). The statistical errors (not shown) are smaller than 3%. 



of the rapidity distribution: 



+ 



0": 



41n(A/sj^j^/2mN) 



erf. 



+ 1 



(45) 



The energy dependence of the sound velocities extracted from the data using Eq. (45) 



is presented in Fig. 15 right. The sound velocities exhibit a clear minimum (usually called 
the softest point) around a beam energy of \/sr 



= 7.6 GeV {30 A GeV). 
As discussed previously the weakening of the transverse and longitudinal expansion 
is expected within the SMES at low SPS energies due to the onset of deconfinement 
which softens the EoS at the early stage. Generally, a softening of the equation of state 
was predicted as a signal for the mixed phase at the transition energy from hadronic to 
partonic matter [HI US, US]- Therefore, we conclude that the data on rapidity spectra 
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of negatively charged pions are indeed compatible with the assumption of the onset of 
deconfinement at the low SPS energies. 



4.5 The Shark Fin and the Tooth 

As discussed in Sec. 3.4 the event-by-event fluctuations of the energy used for particle 
production should lead to fluctuations which are sensitive to the onset of deconfinement. 

The NA49 Collaboration looked for the shark fin structure in the energy dependence 
of the scaled variance of multiplicity distributions in central Pb+Pb collisions [17]. The 
predicted [23] increase of the scaled variance of multiplicity distribution in the NA49 
acceptance by about 0.01 due to the onset of deconfinement is smaller than the systematic 
error on the measurement. Therefore these data can neither support nor disprove the shark 
fin prediction. 



The tooth structure in the energy dependence of Rg/e shown in Fig. 10 right might 
be seen in the event-by-event fluctuations of the K/ir ratio. The energy dependence of 
the fluctuations of this ratio in central Pb+Pb collisions was studied by NA49 using the 
so-called adyn measure [H]. The 'dynamical' K/n fluctuations increase significantly with 
decreasing energy below AOA GeV. It is unclear whether this increase is related to the 
rapid increase of the Rg/e measure predicted due to the onset of deconfinement at energies 
below 30A GeV. 



5 Alternative Approaches 

Several other analyzes of the energy dependence of hadron production properties in central 
Pb+Pb and Au+Au collisions within various theoretical approaches support the hypoth- 
esis that the onset of deconfinement is located at low SPS energies. In particular such 
a result was obtained from studies of hadron yields within a non-equilibrium hadron gas 
model [in] and from an analysis of the time evolution of the relative strangeness yield 
using the momentum integrated Boltzmann equation [SO] . Moreover, a simultaneous anal- 
ysis of the two-pion correlation function and the transverse mass spectra found a plateau 
in the averaged phase-space density at SPS energies which may be associated with the 
onset of deconfinement [5T]. 

Recently a parity violating signal was observed in three-particle azimuthal correlations 
at RHIC [S2]- Such an effect was predicted [S3] since metastable domains may form in 
a QGP where parity and time-reversal symmetries are locally violated. The effect is 
expected to disappear when no QGP is produced in the collisions. It can therefore serve 
as an another indicator for the onset of deconfinement. 

Numerous models have been developed to explain hadron production in reactions of 
heavy nuclei without explicitly invoking a transient QGP phase. The simplest one is the 
statistical hadron gas model [M] which assumes that the hadrochemical freeze-out creates 
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a hadron gas in equilibrium. The temperature, the baryon chemical potential, and the 
hadronization volume are free parameters of the model and are fitted to the data at each 
energy. In this formulation, the hadron gas model cannot make any prediction about the 
energy dependence of hadron production so that an extension of the model was proposed, 
in which the values of the temperature and baryon chemical potential evolve smoothly 
with collision energy [55]. By construction (fit to the energy dependence), the prevailing 
trend in the data is reproduced by the model but important details are not, e.g. the 
decrease of the {K~^) / (ti^) ratio between a/s^y^ = 7.6 and 12.3 GeV {30A and 80A GeV) 
is not well described. The measured ratio of strangeness to pion yield in central Pb+Pb 
collisions at A/s^vAf ~ ^'^■'^ GeV (158^4 GeV) is about 25% lower than the expectation 
for the fully equilibrated hadron gas |55l |56]. Two strategies are followed in order to 
improve the quality of the hadron gas model fits. First, additional parameters have been 
introduced which allow for deviations from equilibrium [STl EHl SH] • Obviously the non- 
equilibrium hadron gas models [S3 ESI ES] with all the parameters fitted separately to the 
data at each energy describe the experimental results significantly better. Secondly the 
equilibrium model was extended to include hypothetical high mass resonance states [59] . 
Again by adding additional free parameters (mass dependence of the resonance state 
density and their branching ratios) the fit quality can be improved. Interestingly, the 
energy dependence of the parameters obtained within the extended hadron gas models is 
interpreted [19],|59] as an indication for the onset of deconfinement at a/s^vtv ~ '^■^ GeV 
{30 A GeV). 

Dynamical models of nucleus-nucleus collisions, such as RQMD [60] , UrQMD [61] and 
HSD [62] treat the initial nucleon-nucleon interactions within a string-hadronic framework. 
In addition these models include effects such as string-string interactions and hadronic re- 
scattering which are expected to be relevant in nucleus-nucleus collisions. RQMD [601163] . 
UrQMD [HH EH |65] and HSD [65] models, like the hadron gas model, fail to describe the 
rapid change of hadron production properties with collision energy in the low SPS energy 
range. 

It was shown that the maximum in relative strangeness production can be reproduced 
by invoking an unusually long lifetime of the fireball at low SPS energies which decreases 
with collision energy [66]. This assumption is however difficult to justify by dynamical 
models of the collision process [611 ES] , and conflicts with the measured energy dependence 
of the two-pion correlation function [67] 168] . 

The step-like structure in the energy dependence of the inverse slope parameter of the 
transverse mass spectra was obtained within the hydrodynamical model by introduction 
of a rapid change of the freeze-out conditions at low SPS energies [6^. However, this 
assumption does not explain the increase of the inverse slope parameter suggested by the 
RHIC results. 

In summary, one is led to conclude that models which do not invoke the onset of 
deconfinement at low SPS energies cannot explain comprehensively and consistently the 
energy dependence of hadron production properties in central Pb+Pb (Au-|-Au) collisions. 
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6 Open problems 



Open questions related to the onset of deconfinement are discussed in this section. First 
the theoretical questions are addressed, then the experimental issues are discussed. 

6.1 Theoretical problems 

Quantum chromodynamics, the commonly accepted theory of strong interactions, was 
developed to model the interactions of quarks and gluons as well as their color neutral 
composites, the hadrons. Thus, in principle, this theory should be able to predict whether 
and via which observables the onset of deconfinement can be experimentally observed in 
nucleus-nucleus collisions. In QCD the strength of the strong force between two quarks 
increases with their distance. This property of the theory has as a consequence that, in 
general, predictions are either very difficult or impossible to calculate. Presently, there are 
no quantitative QCD predictions yet concerning the onset of deconfinement in nucleus- 
nucleus collisions. 

The SMES model |1] assumes statistical particle production at the early stage of 
nucleus-nucleus collisions. The data on nucleus-nucleus collisions at RHIC, in particular 
results on anisotropic fiow, seem to require a large degree of equilibration at the early 
stage of collisions [70]. Thus the assumption of statistical particle production received 
independent experimental support. 

The additional assumptions of the SMES, which lead to the kink and horn predictions, 
concern entropy and strangeness conservation during expansion and freeze-out. In fact, 
the model predictions remain at least qualitatively unchanged if one assumes an approxi- 
mate proportionality of the final state entropy and strangeness to their early stage values. 
There is no easy way to prove or disprove these weaker requirements. 

The predictions concerning the step and dale require assumptions concerning the equa- 
tion of state of strongly interacting matter. In particular, they rely on a general feature 
of the EoS, i.e. the existence of the softest point when the transition between QGP and 
confined matter occurs. Thus, the qualitative predictions are independent of the nature of 
the transition (e.g. cross-over or 1** order phase transition, full or partial energy stopping 
and thermalization) . However, the quantitative predictions are sensitive and a consistent 
description of the hydrodynamical evolution has not been achieved yet. In particular, the 
bag model equation of state with a 1'^* order phase transition to hadron gas leads to a 
significant over-prediction of transverse flow [71] . Further studies are needed. 

The predictions concerning the tooth and shark fin are derived assuming that the 
early stage volume fluctuations can be neglected for central collisions. It is unclear, to 
which extent this condition is consistent with recent attempts to describe multiplicity 
distributions and high transverse momentum spectra in p+p interactions by a statistical 
model with volume fluctuations [72] . Further studies are needed. 
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6.2 Experimental issues 

The experimental results indicating the onset of deconfinement were obtained mainly by 
the NA49 experiment at the CERN SPS. Clearly, a confirmation of these measurements 
is necessary. 

A beam energy scan program at the BNL RHIC has recently started. Pilot results at 
y^NN ~ 9-2 and 20 GeV are in agreement with the corresponding NA49 data [73]. New 
RHIC data being taken by the STAR experiment in 2010 will allow a more conclusive 
verification of the NA49 results. 

At the CERN SPS the new experiment NA61 started a two dimensional system size and 
beam energy scan in 2009, which will continue over the next 4 years. The measurements 
aim to verify the existence of the onset of deconfinement in collisions of medium size 
nuclei. Moreover, they will allow to study the expected disappearance of the signals in 
collisions of light nuclei. 

7 Summary and Conclusions 

In this review we present the experimental and theoretical status of the evidence for the 
threshold of quark-gluon plasma creation in high energy nucleus-nucleus interactions. The 
location in energy of this so-called onset of deconfinement, as well as key experimental 
signals were predicted by the statistical model of the early stage of the collision process [1] . 
These signals were searched for and observed within the energy scan program of the NA49 
Collaboration at the CERN SPS. Together with measurements at lower (LBL, JINR, SIS, 
BNL AGS) and higher (BNL RHIC) energies the properties of hadron production in heavy 
ion collisions were established in a broad energy range. Their energy dependence led to the 
conclusion that the predicted signals of the onset of deconfinement appear in a common 
energy domain covered by the SPS at CERN. These features of the data serve as strong 
experimental evidence for the existence of the onset of deconfinement and thus for the 
quark-gluon plasma itself. 

Quantitative model predictions, discussed in this review, are derived within the statis- 
tical approach to particle production in high energy collisions. The use of this approach 
has a two-fold justification. First, it naturally includes the concept of phases of strongly 
interacting matter and the transition between them. Second, it is successful in describing 
numerous features of the experimental data. The relation between the statistical approach 
and the commonly accepted theory of strong interactions, QCD, remains unclear. This 
is because QCD is difficult or impossible to evaluate in the energy region relevant for 
multi-particle production in general and for the phase transitions of strongly interacting 
matter in particular. 

New experimental programs have started at the CERN SPS and BNL RHIC which are 
devoted to the study of nucleus-nucleus collisions in the energy region where the NA49 
experiment found evidence for the onset of deconfinement. The STAR experiment at 
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RHIC will provide a necessary confirmation of these results. The new CERN experiment 
NA61 will address the questions how this observed phenomenon depends on the volume 
of matter and what the properties of the transition region are. 
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8 Appendices 



8.1 Strangeness Enhancement and J/ifj Suppression 

The idea of strangeness enhancement as a quark-gluon plasma signal in nucleus-nucleus 
(A+A) collisions was formulated a long time ago [8]. It was based on the estimate that the 
strangeness equilibration time in the QGP is of the same order (ca 10 fm/c) as the expected 
life time of the fireball formed in A+A collisions. Thus in the case of QGP creation 
strangeness is expected to approach its QGP equilibrium value. This equilibrium value 
is significantly higher than the level of strangeness production in nucleon-nucleon (N+N) 
collisions. Strangeness production in secondary hadronic interactions was estimated to be 
negligibly small. Therefore, if QGP is not formed, strangeness yields would be expected to 
be much lower than those predicted by equilibrium QGP calculations. Thus at that time 
a simple and elegant signature of QGP creation appeared: a transition to QGP should be 
signaled by an increase of the strangeness production level to the QGP equilibrium value. 

The actual study has been done in the following way. The strangeness to pion ratio 
quantified by the ratios, 

(A) + {K + K) {K+) 
(tt) (7r+) 

was measured and analyzed. One expected that the ratios should increase strongly in A+A 
collisions if the QGP was formed. To reveal the specific increase of the strangeness/pion 
ratio in A+A collisions due to QGP formation the strangeness enhancement factor was 
introduced: 

where superscripts and correspond respectively to A+A and N+N collisions at the 
same N+N cm. energy y/sj^^. The confrontation of this expectation with data was for the 
first time possible in 1988 when results from ^^S and ^^Si beams at the SPS and the AGS 
became available. Experiment NA35 reported that in central S+S collisions at 200^4 GeV 
the strangeness to pion ratio is indeed about 2 times higher than in N+N interactions at 
the same energy per nucleon [71]. But an even larger enhancement {Rs is about 5) was 
measured by E802 in Si+A collisions at the AGS. The data on central Au+Au collisions 
at low AGS energies 2y4-10A GeV completed the picture: strangeness enhancement is 



observed at all energies, and it is stronger at lower energies, i.e. the ratio Rg Eq.(47) 



increases monotonically with decreasing \/s^j^. Figure 16 shows a compilation of recent 
data ^75j . 



Moreover, the enhancement factor Rg Eq.(47) should evidently go to infinity at the 
threshold energy of strange hadron production in N+N collisions. At low AGS energies 
one does not expect the creation of a QGP and therefore the substantial strangeness 
enhancement is clearly of a different origin. The low level of strangeness production 



38 



15 




Figure 16: Energy dependence of strangeness enhancement in central Pb+Pb (Au+Au) 
collisions [75] . 



in p+p interactions as compared to the strangeness yield in central A+A collisions can 
be understood to a large extent in the statistical model as due to the effect of exact 
strangeness conservation. The canonical ensemble treatment of strangeness conservation 
leads to additional suppression factors imposed on strange hadron production in small 
systems such as created in p+p collisions. In any case, the AGS measurements indicating 
a strangeness enhancement Rs larger than that at the SPS show clearly that the simple 
concept of strangeness enhancement as a signal of QGP production does not work. 

The standard picture of J/ip production in collisions of hadrons and nuclei assumes 
a two step process: the creation of a cc pair in hard parton collisions at the very early 
stage of the reaction and subsequent formation of a bound charmonium state. It was 
proposed [U] to use the J/?/' as a probe for deconfinement in the study of A+A collisions. 
The argument was that in a QGP color screening dissolves initially created J/ip mesons 
into c and c quarks which at hadronization form open charm hadrons. As the initial yield 
of J/tp is believed to have the same A-dependence as the yield of Dell-Yan lepton pairs, 
the measurement of a weaker A-dependence of the final J/ip yield {J/ip suppression) 
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would signal charmonium absorption and therefore creation of QGP. 

Production of charmonium states J/ip and ip' was measured in A+A collisions at the 
CERN SPS over a period of 15 years by the NA38, NA50 and NA60 Collaborations [IS]. 
The A-dependence of J /ip production in p+A is weaker than A^ (approximately A°'^). It 
was suggested that this Jjip suppression is due to absorption on nucleons in the target 
nucleus. The data on oxygen and sulfur collisions on nuclei at 200A GeV also indicated 
considerable suppression. To improve the fit of the data a new source of J /ip absorption 
was introduced: the absorption on hadronic secondaries ('co-movers'). Finally in central 
Pb+Pb collisions at 158A GeV the measured suppression was found to be significantly 
stronger than expected in the models including nuclear and co-mover suppression. This 
'anomalous' J/ip suppression was interpreted as evidence of QGP creation in Pb+Pb 
collisions at the CERN SPS. The uncertainties in estimates of the J /if) absorption by 
target nucleons and co- movers make J/ip suppression a problematic QGP signal. An 
essential part of the J /ip yield comes from decays of excited charmonium states like ip' 
and X- All of them have different melting temperatures and absorption cross-sections. 

Alternative approaches have been developed, namely the statistical [77] and the sta- 
tistical coalescence [751 [79] models of J /ip production, which reproduce the A-dependence 
of the J/ip yield at SPS energies reasonably well. They are based on different physics 
pictures than the one leading to J/ip suppression as the signal of quark-gluon plasma 
creation. Specifically, the statistical model [77] assumes statistical production of J/ip 
mesons at hadronization, whereas in the coalescence model statistical coalescence of c 
and c quarks at hadronization is assumed [751 [79] . In both models the J /ip yield is 
neither related to the J/psi suppression in the hadron gas nor in the quark-gluon plasma. 
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8.2 Main carriers of s and s quarks 



main strangeness carriers 
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Figure 17: Main carriers of s and s quarks. K~ and as well as A and A yields are 
sensitive to the strangeness content and baryon density. and yields are sensitive 
mainly to strangeness content. 



The distribution of s and s quarks among the most abundantly produced hadrons is 
considered here for the case of nucleus-nucleus collisions at the SPS energies. The sketch 
presented in Fig. 17 illustrates the following discussion. 

The colliding nuclei have net numbers of s and s quarks equal to zero. As strangeness 
is conserved in strong interactions the numbers of s and s quarks in the final state have 
to be equal. 

Kaons are the lightest strange hadrons. The isospin partners and K'^ mesons carry 
s quarks, whereas K~ and carry s quarks. The (approximate) symmetry of the initial 
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state and isospin conservation in strong interactions imply that: 



{K+) ^ 



(48) 



and 



(49) 



The s-quarks are also carried by the hghtest anti-baryon, A. Its fraction is however small 
(less than 5%) at the AGS and SPS energies due to suppression of the anti-baryon yield 
by the high net-baryon density. Consequently, and mesons carry each about half 
of all the anti-strange quarks produced in A-l-A collisions at AGS and SPS energies. Thus, 
their yields are nearly proportional to the total number of produced s and s quarks. 

This is not the case for K~ and mesons. A significant fraction of s-quarks (about 
50% in central Pb+Pb collisions at 158 A-GeV) is carried by hyperons. In addition 
this fraction strongly depends on collision energy. Consequently, the fraction of s-quarks 
carried by anti-kaons, and K , is also dependent on collision energy and cannot be 
used easily to quantify strangeness production. In the Es ratio all main carriers of strange 
and anti-strange quarks are included. The neglected contribution of A and other hyperons 
and anti-hyperons is about 10% at SPS energies. Both the {K^) / (7r+) and Es ratios are 
approximately, within 5% at SPS energies, proportional to the ratio of total multiplicity 
of s and s quarks to the multiplicity of pions. It should be noted that the {K^) / {ji^) 
ratio is expected to be similar (within about 10%) for n + p, and n + n interactions 

at 158^4 GeV, whereas the Eg ratio is independent of the isospin of nucleon-nucleon 
interactions. 
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8.3 Onset of deconfinement and critical point 




Figure 18: Onset of deconfinement and critical point. 



This appendix discusses relations between the onset of deconfinement, the critical 
point of strongly interacting matter and the possibilities of their experimental study in 



relativistic nucleus-nucleus collisions. The two sketches presented in Fig. 18 should help 
to understand the basic ideas. 

The onset of deconfinement refers to the beginning of the creation of a deconfined 
state of strongly interacting matter (ultimately a quark-gluon plasma) at the early stage 
of nucleus-nucleus collisions when increasing the collision energy. With increasing collision 
energy the energy density of matter created at the early stage of A+A collisions increases. 
Thus, if there are two phase^ of matter separated by the transition region (solid and 

^ The discussed two phase diagram is the simplest one which aUows to introduce the concepts of the 
onset of deconfinement and the transition region. There are numerous suggestions of phase diagrams 
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dotted lines) as indicated in Fig. 18 left the early stage (the red point) first has to hit 
and then move above the transition region. Therefore, the existence of the onset of 
deconfinement is the most straightforward consequence of the existence of two phases of 
strongly interacting matter, i.e. confined matter and QGP. The experimental observation 
of the onset of deconfinement required a one dimensional scan in collision energy with 
heavy ions as performed by NA49. All signals of the onset of deconfinement discussed in 
this paper relate to the difference in properties of confined matter and QGP. They are 
not sensitive to the structure of the transition region. 

Discovery of the onset of deconfinement implies the existence of QGP and of a tran- 
sition region between confined and QGP phases. Numerous possibilities concerning the 
structure of the transition region are under discussion (see e.g., Ref. [81]). The most 



popular one |_82j, sketched in Fig. 18, claims that a 1*** order phase transition (thick gray 
line) separates both phases in the high baryonic chemical potential domain. In the low 
baryonic chemical potential domain a rapid crossover is expected (dotted line). The end 
point of the 1'^* order phase transition line is the critical point. 

The characteristic signatures of the critical point can be observed if the freeze-out 
point (blue square in Fig. 18 right) is located close to the critical point. The analysis of 
the existing experimental data |5H] indicates that the location of the freeze-out point in 
the phase diagram depends on the collision energy and the mass of the colliding nuclei. 
This dependence is schematically indicated in Fig. 18 right. Thus the experimental search 
for the critical point requires a two-dimensional scan in collision energy and size of the 
colliding nuclei. The NA61 experiment [T71 118] at the CERN SPS started this scan in 
2009. It should be completed within several years. Note, that a two dimensional scan is 
actually required for any study of the structure of the transition region, independent of 
the hypothesis tested. 

The transition region can be studied experimentally in nucleus-nucleus collisions only 
at T,fiB values which correspond to collision energies higher than the energy of the onset of 



deconfinement. This important conclusion is easy to understand when looking at Fig 18 



Signals of the critical point can be observed provided the freeze-out point is close to it 
(see Fig. 18 right). On the other hand, by definition the critical point is located on the 
transition line. Furthermore, the energy density at the early stage of the collision is, of 
course, higher than the energy density at freeze-out. Thus, the condition that the freeze- 
out point is near the critical point implies that the early stage of the system is above (or 
on) the transition line. This in turn means that the optimal energy range for the search 
for the critical point lies above the energy of the onset of deconfinement (see Fig. 18 left). 
This general condition limits the search for the critical point to the collision energy range 
Elab > 30A GeV. 



with a significantly richer structure (see e.g., Ref. [8U]). 
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